This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.
Issue Howard Bayes
Time series 1999 - present 1977 - present
Bias in SWHS Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests. Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data. Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas. Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with. By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.
**Figure 1.**- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.


**Figure 2.**- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).


Results

**Figure X.**- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

**Figure 3.**- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 3.**- Total rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

**Figure 4.**- Black rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


And black rockfish releases…

**Figure 5.**- Black rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 6.**- Yellow rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 7.**- Yellow rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 8.**- DSR rockfish (excluding yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 9.**- DSR rockfish releases (including yelloweye) 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 11.**- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 12.**- Slope rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

**Figure 13.**- *Other* rockfish (non-pelagic and non-yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 14.**- *Other* rockfish (non-pelagic and non-yelloweye) releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



Total Biomass Removal Estimates

**Figure 15.**- Black rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.



**Figure 16.**- Yellow rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

**Figure 17.**- Pelagic rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 18.**- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 19.**- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


Model fit

Logbook residuals

**Figure 20.**- Residuals from logbook harvests.

Figure 20.- Residuals from logbook harvests.


SWHS residuals

**Figure 21.**- Residuals from SWHS harvests.

Figure 21.- Residuals from SWHS harvests.



**Figure 22.**- Residual of SWHS releases.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

**Figure 23.**- Mean percent of harvest by charter anglers.

Figure 23.- Mean percent of harvest by charter anglers.


When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

**Figure 24.**- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.


**Figure 25.**- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 26.**- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 27.**- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.


**Figure 28.**- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.


**Figure 29.**- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

**Figure 30.**- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.


Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

**Figure 31.**- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

**Figure 32.**- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

**Figure 33.**- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

**Figure 34.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

**Figure 35.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

**Figure 36.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.



P(slope|non-pelagic & non-yellowye) For release estimates

## NULL



Weight Fits

**Figure 37.**- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 38.**- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 39.**- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 40.**- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 41.**- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.5046296

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.15 while 0.85 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable

area_n

year

user

mean

median

lower

upper

Rhat

sd

skew

area

p_pelagic

11

1

1

0.467

0.467

0.228

0.712

1.011

0.126

0.000

CSEO

p_pelagic

11

11

2

0.404

0.399

0.186

0.652

1.015

0.121

0.041

CSEO

p_pelagic

11

15

1

0.430

0.426

0.198

0.684

1.017

0.124

0.032

CSEO

p_pelagic

11

15

2

0.394

0.389

0.172

0.646

1.015

0.122

0.041

CSEO

p_pelagic

11

17

2

0.401

0.395

0.182

0.640

1.014

0.120

0.050

CSEO

p_pelagic

11

2

2

0.400

0.397

0.183

0.652

1.013

0.121

0.025

CSEO

p_pelagic

11

22

2

0.410

0.403

0.192

0.659

1.013

0.120

0.058

CSEO

p_pelagic

11

4

1

0.468

0.466

0.232

0.709

1.012

0.126

0.016

CSEO

p_pelagic

11

5

2

0.399

0.392

0.183

0.656

1.011

0.123

0.057

CSEO

p_pelagic

11

7

1

0.468

0.466

0.236

0.713

1.012

0.125

0.016

CSEO

p_pelagic

11

8

1

0.467

0.466

0.228

0.712

1.011

0.127

0.008

CSEO

p_pelagic

11

9

1

0.467

0.467

0.230

0.713

1.017

0.125

0.000

CSEO

p_pelagic

13

1

1

0.301

0.281

0.101

0.595

1.014

0.132

0.152

NSEI

p_pelagic

13

1

2

0.123

0.107

0.035

0.313

1.014

0.071

0.224

NSEI

p_pelagic

13

10

2

0.131

0.115

0.037

0.308

1.016

0.073

0.220

NSEI

p_pelagic

13

11

2

0.131

0.114

0.039

0.320

1.019

0.074

0.230

NSEI

p_pelagic

13

13

2

0.137

0.121

0.041

0.315

1.013

0.073

0.219

NSEI

p_pelagic

13

14

1

0.310

0.298

0.119

0.566

1.010

0.119

0.101

NSEI

p_pelagic

13

15

1

0.284

0.268

0.106

0.554

1.010

0.117

0.137

NSEI

p_pelagic

13

15

2

0.131

0.117

0.041

0.295

1.010

0.067

0.208

NSEI

p_pelagic

13

2

1

0.300

0.281

0.098

0.585

1.015

0.131

0.145

NSEI

p_pelagic

13

2

2

0.124

0.105

0.033

0.312

1.014

0.073

0.259

NSEI

p_pelagic

13

3

1

0.301

0.282

0.100

0.598

1.015

0.132

0.144

NSEI

p_pelagic

13

4

1

0.302

0.283

0.108

0.588

1.010

0.128

0.148

NSEI

p_pelagic

13

4

2

0.125

0.110

0.033

0.310

1.018

0.073

0.204

NSEI

p_pelagic

13

5

2

0.124

0.107

0.034

0.311

1.012

0.072

0.236

NSEI

p_pelagic

13

6

2

0.126

0.109

0.034

0.315

1.012

0.073

0.231

NSEI

p_pelagic

13

7

2

0.126

0.111

0.035

0.308

1.015

0.072

0.210

NSEI

p_pelagic

13

8

1

0.308

0.292

0.108

0.598

1.012

0.128

0.125

NSEI

p_pelagic

13

9

2

0.129

0.111

0.036

0.308

1.015

0.073

0.248

NSEI

p_pelagic

16

1

1

0.400

0.392

0.158

0.673

1.032

0.138

0.058

SSEO

p_pelagic

16

1

2

0.376

0.369

0.142

0.654

1.033

0.136

0.052

SSEO

p_pelagic

16

10

1

0.403

0.397

0.160

0.676

1.040

0.137

0.044

SSEO

p_pelagic

16

10

2

0.382

0.372

0.149

0.657

1.033

0.135

0.074

SSEO

p_pelagic

16

11

1

0.404

0.400

0.163

0.685

1.041

0.136

0.029

SSEO

p_pelagic

16

11

2

0.385

0.379

0.147

0.660

1.034

0.134

0.045

SSEO

p_pelagic

16

12

1

0.401

0.397

0.159

0.670

1.028

0.135

0.030

SSEO

p_pelagic

16

12

2

0.384

0.377

0.152

0.661

1.022

0.133

0.053

SSEO

p_pelagic

16

13

1

0.405

0.398

0.162

0.674

1.030

0.135

0.052

SSEO

p_pelagic

16

13

2

0.383

0.376

0.148

0.646

1.024

0.132

0.053

SSEO

p_pelagic

16

14

1

0.391

0.384

0.156

0.655

1.027

0.130

0.054

SSEO

p_pelagic

16

14

2

0.382

0.373

0.160

0.648

1.030

0.130

0.069

SSEO

p_pelagic

16

15

1

0.396

0.391

0.163

0.657

1.034

0.131

0.038

SSEO

p_pelagic

16

15

2

0.386

0.376

0.156

0.656

1.027

0.133

0.075

SSEO

p_pelagic

16

16

1

0.391

0.382

0.163

0.659

1.022

0.130

0.069

SSEO

p_pelagic

16

16

2

0.384

0.377

0.155

0.654

1.017

0.130

0.054

SSEO

p_pelagic

16

17

1

0.446

0.441

0.208

0.710

1.020

0.131

0.038

SSEO

p_pelagic

16

17

2

0.399

0.392

0.171

0.663

1.037

0.131

0.053

SSEO

p_pelagic

16

18

1

0.426

0.419

0.192

0.684

1.020

0.129

0.054

SSEO

p_pelagic

16

18

2

0.399

0.395

0.170

0.663

1.022

0.128

0.031

SSEO

p_pelagic

16

19

1

0.420

0.414

0.188

0.687

1.015

0.130

0.046

SSEO

p_pelagic

16

19

2

0.399

0.394

0.170

0.659

1.022

0.128

0.039

SSEO

p_pelagic

16

2

1

0.400

0.394

0.148

0.682

1.028

0.138

0.043

SSEO

p_pelagic

16

2

2

0.378

0.372

0.137

0.644

1.035

0.135

0.045

SSEO

p_pelagic

16

20

1

0.414

0.410

0.190

0.665

1.019

0.124

0.032

SSEO

p_pelagic

16

20

2

0.404

0.398

0.175

0.671

1.013

0.128

0.047

SSEO

p_pelagic

16

21

1

0.398

0.389

0.180

0.659

1.014

0.124

0.073

SSEO

p_pelagic

16

21

2

0.399

0.393

0.184

0.655

1.016

0.125

0.048

SSEO

p_pelagic

16

22

2

0.405

0.399

0.187

0.664

1.014

0.124

0.048

SSEO

p_pelagic

16

3

1

0.398

0.395

0.149

0.676

1.037

0.139

0.022

SSEO

p_pelagic

16

3

2

0.376

0.369

0.137

0.650

1.037

0.136

0.051

SSEO

p_dsr

16

31

1

0.407

0.407

0.376

0.437

1.016

0.016

0.000

SSEO

p_pelagic

16

4

1

0.404

0.403

0.155

0.677

1.031

0.141

0.007

SSEO

p_pelagic

16

4

2

0.379

0.378

0.136

0.648

1.029

0.136

0.007

SSEO

p_pelagic

16

5

1

0.402

0.397

0.148

0.670

1.040

0.138

0.036

SSEO

p_pelagic

16

5

2

0.378

0.371

0.140

0.659

1.034

0.136

0.051

SSEO

p_pelagic

16

6

1

0.399

0.396

0.149

0.678

1.030

0.139

0.022

SSEO

p_pelagic

16

6

2

0.380

0.371

0.142

0.654

1.029

0.135

0.067

SSEO

p_pelagic

16

7

1

0.402

0.399

0.152

0.677

1.045

0.139

0.022

SSEO

p_pelagic

16

7

2

0.376

0.370

0.137

0.656

1.034

0.135

0.045

SSEO

p_pelagic

16

8

1

0.405

0.400

0.154

0.670

1.043

0.137

0.036

SSEO

p_pelagic

16

8

2

0.380

0.373

0.148

0.655

1.038

0.135

0.052

SSEO

p_pelagic

16

9

1

0.403

0.399

0.155

0.682

1.038

0.140

0.029

SSEO

p_pelagic

16

9

2

0.383

0.380

0.137

0.660

1.036

0.138

0.022

SSEO

p_pelagic

4

1

1

0.537

0.548

0.261

0.788

1.010

0.137

-0.080

PWSO

p_pelagic

4

1

2

0.544

0.553

0.271

0.795

1.012

0.138

-0.065

PWSO

p_pelagic

4

10

1

0.536

0.541

0.271

0.780

1.013

0.135

-0.037

PWSO

p_yellow

4

10

1

0.878

0.893

0.694

0.975

1.011

0.075

-0.199

PWSO

p_pelagic

4

11

1

0.536

0.541

0.266

0.780

1.011

0.136

-0.037

PWSO

p_yellow

4

11

2

0.875

0.890

0.693

0.975

1.011

0.073

-0.204

PWSO

p_pelagic

4

12

2

0.545

0.552

0.267

0.791

1.012

0.137

-0.051

PWSO

p_yellow

4

12

1

0.878

0.893

0.689

0.972

1.016

0.072

-0.208

PWSO

p_pelagic

4

14

1

0.531

0.536

0.267

0.772

1.012

0.131

-0.038

PWSO

p_pelagic

4

14

2

0.543

0.548

0.265

0.796

1.014

0.135

-0.037

PWSO

p_yellow

4

14

1

0.878

0.892

0.702

0.975

1.024

0.074

-0.190

PWSO

p_pelagic

4

15

2

0.545

0.549

0.269

0.782

1.014

0.132

-0.030

PWSO

p_yellow

4

15

1

0.878

0.896

0.690

0.971

1.012

0.073

-0.246

PWSO

p_pelagic

4

16

1

0.534

0.540

0.268

0.789

1.010

0.136

-0.044

PWSO

p_yellow

4

16

2

0.875

0.892

0.685

0.971

1.013

0.074

-0.230

PWSO

p_pelagic

4

17

2

0.549

0.556

0.279

0.782

1.010

0.133

-0.053

PWSO

p_pelagic

4

19

1

0.537

0.538

0.275

0.781

1.013

0.131

-0.008

PWSO

p_pelagic

4

2

1

0.535

0.540

0.249

0.781

1.016

0.136

-0.037

PWSO

p_pelagic

4

2

2

0.542

0.547

0.273

0.790

1.014

0.134

-0.037

PWSO

p_pelagic

4

3

2

0.542

0.548

0.260

0.792

1.011

0.138

-0.043

PWSO

p_yellow

4

3

1

0.880

0.895

0.696

0.979

1.012

0.073

-0.205

PWSO

p_pelagic

4

4

1

0.537

0.539

0.265

0.776

1.014

0.136

-0.015

PWSO

p_yellow

4

4

1

0.881

0.897

0.700

0.979

1.010

0.073

-0.218

PWSO

p_pelagic

4

5

1

0.535

0.541

0.261

0.786

1.011

0.135

-0.044

PWSO

p_pelagic

4

5

2

0.541

0.545

0.255

0.791

1.014

0.138

-0.029

PWSO

p_pelagic

4

6

1

0.536

0.541

0.261

0.784

1.013

0.137

-0.037

PWSO

p_pelagic

4

6

2

0.545

0.552

0.258

0.791

1.011

0.136

-0.051

PWSO

p_pelagic

4

7

1

0.537

0.540

0.265

0.786

1.018

0.137

-0.022

PWSO

p_pelagic

4

7

2

0.544

0.550

0.265

0.783

1.011

0.135

-0.045

PWSO

p_pelagic

4

8

1

0.540

0.542

0.265

0.792

1.010

0.137

-0.015

PWSO

p_yellow

4

8

2

0.876

0.890

0.686

0.975

1.011

0.075

-0.186

PWSO

p_pelagic

4

9

2

0.544

0.549

0.283

0.785

1.013

0.134

-0.037

PWSO

p_black

8

40

1

0.461

0.414

0.340

0.957

1.021

0.160

0.293

afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable

area_n

year/species

mean

median

lower

upper

Rhat

sd

skew

year

area

beta0_pH

11

4

0.783

0.914

-1.487

1.902

1.268

0.764

-0.171

CSEO

beta1_pH

11

4

2.010

1.918

0.715

4.221

1.213

0.792

0.116

CSEO

beta0_pelagic

11

-0.137

-0.111

-0.604

0.148

1.130

0.188

-0.138

CSEO

beta0_pH

12

4

-0.146

0.914

-17.076

2.226

1.472

4.209

-0.252

EWYKT

beta1_pH

12

4

4.233

3.160

1.802

20.983

1.469

4.216

0.254

EWYKT

beta3_pH

12

4

41.902

42.116

40.206

43.764

1.156

0.939

-0.228

EWYKT

beta0_pH

13

4

0.401

0.736

-3.051

1.586

1.374

1.116

-0.300

NSEI

beta1_pH

13

4

1.593

1.290

0.313

5.004

1.383

1.118

0.271

NSEI

beta0_pH

14

4

0.310

0.435

-1.491

1.276

1.375

0.698

-0.179

NSEO

beta1_pH

14

4

2.056

1.942

0.905

3.976

1.363

0.769

0.148

NSEO

beta0_yellow

14

-0.298

-0.236

-1.341

-0.010

1.268

0.350

-0.177

NSEO

beta0_pH

15

4

0.486

0.684

-1.237

1.532

1.364

0.732

-0.270

SSEI

beta1_pH

15

4

1.515

1.359

0.323

3.261

1.306

0.765

0.204

SSEI

beta3_pH

15

4

41.309

41.266

38.831

43.706

1.203

1.200

0.036

SSEI

beta0_pH

16

4

0.761

0.907

-1.699

2.152

1.344

0.963

-0.152

SSEO

beta1_pH

16

4

3.619

3.460

2.104

6.231

1.324

1.013

0.157

SSEO

beta0_pelagic

16

-0.452

-0.361

-1.340

0.071

1.133

0.384

-0.237

SSEO

mu_beta2_pH

2

1

1.449

0.766

0.073

7.081

1.131

2.083

0.328

NG

mu_beta0_pH

3

4

0.442

0.773

-3.752

1.657

1.433

1.356

-0.244

PWSI

mu_beta1_pH

3

4

2.274

2.031

0.642

5.233

1.230

1.312

0.185

PWSI

beta2_pH

6

1

1.312

0.468

0.108

8.151

1.122

2.231

0.378

SOKO2SAP

beta2_pH

8

1

1.108

0.390

0.114

7.392

1.116

1.841

0.390

afognak

beta3_pH

9

1

31.955

31.625

27.153

39.698

1.114

2.817

0.117

eastside

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.01

307.00

193.00

944.00

332.00

Prop. of total

0.61

0.39

0.74

0.26

Prop. median less than 1

0.32

0.23

0.35

0.45

Prop. median less than 10

0.46

0.57

0.58

0.80

Prop. median less than 50

0.67

0.78

0.76

0.86

Median of the median values

15.18

7.21

4.29

1.51

Mean of the median values

234.96

87.28

157.68

90.52

Mean skew of posteriors

0.26

0.30

0.27

0.24

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.1

3.00

4.00

43.00

40.00

Prop. of total

0.43

0.57

0.52

0.48

Prop. median less than 1

1.00

1.00

0.79

0.52

Prop. median less than 10

1.00

1.00

0.86

0.92

Prop. median less than 50

1.00

1.00

0.95

0.95

Median of the median values

0.00

0.54

0.01

0.51

Mean of the median values

0.01

0.53

10.50

10.71

Mean skew of posteriors

0.15

0.20

0.14

0.17

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable

Freq

beta0_pelagic

2

beta0_pH

6

beta0_yellow

1

beta1_pH

6

beta2_pH

2

beta3_pH

3

eps_pH

3

Ho_ay

1

Ho_ayu

4

Hs_ayu

1

Hy_ay

3

Hy_ayu

3

mu_beta0_pH

1

mu_beta1_pH

1

mu_beta2_pH

1

pH

2

R_ayg

1

R_ayu

1

Rb_ayg

1

Rb_ayu

1

Rd_ayg

1

Rd_ayu

1

re_pH

5

Ro_ay

15

Ro_ayg

16

Ro_ayu

28

Rp_ayg

1

Rp_ayu

1

Rs_ay

2

Rs_ayg

8

Rs_ayu

19

Ry_ay

12

Ry_ayg

3

Ry_ayu

24

Parameter estimates:

Summary Table of Parameter Estimates
Parameter mean sd Lower_CI Median Upper_CI
mu_bc_H[1] -0.016 0.069 -0.142 -0.021 0.127
mu_bc_H[2] 0.072 0.055 -0.029 0.071 0.188
mu_bc_H[3] -0.267 0.067 -0.398 -0.268 -0.135
mu_bc_H[4] -0.924 0.180 -1.287 -0.918 -0.584
mu_bc_H[5] 1.313 1.102 0.120 1.063 4.002
mu_bc_H[6] -2.400 0.313 -3.000 -2.413 -1.762
mu_bc_H[7] -0.457 0.103 -0.662 -0.453 -0.266
mu_bc_H[8] 0.272 0.346 -0.308 0.241 1.067
mu_bc_H[9] -0.247 0.125 -0.496 -0.248 0.005
mu_bc_H[10] 0.058 0.066 -0.065 0.057 0.191
mu_bc_H[11] -0.126 0.037 -0.200 -0.125 -0.054
mu_bc_H[12] -0.250 0.107 -0.487 -0.243 -0.053
mu_bc_H[13] -0.160 0.072 -0.297 -0.162 -0.016
mu_bc_H[14] -0.305 0.089 -0.488 -0.303 -0.127
mu_bc_H[15] -0.339 0.047 -0.435 -0.339 -0.246
mu_bc_H[16] -0.141 0.366 -0.785 -0.156 0.643
mu_bc_R[1] 1.042 0.121 0.803 1.042 1.276
mu_bc_R[2] 1.010 0.080 0.842 1.014 1.156
mu_bc_R[3] 0.224 0.090 0.035 0.225 0.399
mu_bc_R[4] 0.678 0.198 0.270 0.686 1.040
mu_bc_R[5] 1.674 0.375 0.933 1.673 2.392
mu_bc_R[6] -1.429 0.339 -2.109 -1.426 -0.737
mu_bc_R[7] 0.509 0.171 0.157 0.510 0.836
mu_bc_R[8] 0.651 0.145 0.385 0.647 0.943
mu_bc_R[9] 0.383 0.167 0.036 0.390 0.694
mu_bc_R[10] 0.897 0.104 0.694 0.899 1.102
mu_bc_R[11] 1.113 0.071 0.981 1.111 1.255
mu_bc_R[12] 0.873 0.162 0.540 0.875 1.176
mu_bc_R[13] 1.031 0.089 0.847 1.031 1.199
mu_bc_R[14] 0.861 0.123 0.620 0.862 1.095
mu_bc_R[15] 0.801 0.097 0.603 0.802 0.987
mu_bc_R[16] 1.353 0.117 1.123 1.356 1.574
tau_pH[1] 2.425 0.183 2.091 2.421 2.806
tau_pH[2] 1.606 0.154 1.323 1.600 1.923
tau_pH[3] 1.196 0.113 0.989 1.191 1.429
tau_pH[4] 1.017 0.141 0.771 1.009 1.318
tau_pH[5] 0.427 0.076 0.295 0.422 0.594
beta0_pH[1,1] 0.549 0.183 0.175 0.553 0.917
beta0_pH[2,1] 1.460 0.176 1.102 1.468 1.772
beta0_pH[3,1] 1.389 0.226 0.860 1.411 1.769
beta0_pH[4,1] 1.601 0.230 1.073 1.615 1.997
beta0_pH[5,1] -0.776 0.414 -1.829 -0.728 -0.120
beta0_pH[6,1] -0.523 0.495 -1.875 -0.377 0.090
beta0_pH[7,1] 0.757 0.011 0.736 0.757 0.778
beta0_pH[8,1] -0.680 0.353 -1.636 -0.629 -0.120
beta0_pH[9,1] -0.351 0.355 -1.122 -0.299 0.179
beta0_pH[10,1] 0.295 0.209 -0.128 0.310 0.683
beta0_pH[11,1] 0.069 0.311 -0.622 0.087 0.623
beta0_pH[12,1] 0.629 0.190 0.228 0.640 0.971
beta0_pH[13,1] -0.654 0.151 -0.970 -0.650 -0.368
beta0_pH[14,1] 0.013 0.181 -0.350 0.016 0.361
beta0_pH[15,1] 0.066 0.193 -0.376 0.087 0.385
beta0_pH[16,1] -0.879 0.221 -1.344 -0.865 -0.477
beta0_pH[1,2] 2.562 0.085 2.396 2.562 2.737
beta0_pH[2,2] 3.079 0.022 3.037 3.079 3.123
beta0_pH[3,2] 3.625 0.041 3.547 3.625 3.705
beta0_pH[4,2] 2.937 0.019 2.899 2.937 2.975
beta0_pH[5,2] 4.550 1.393 2.737 4.268 8.261
beta0_pH[6,2] 3.327 0.152 3.037 3.323 3.633
beta0_pH[7,2] 2.033 0.049 1.939 2.033 2.131
beta0_pH[8,2] 2.823 0.064 2.698 2.823 2.951
beta0_pH[9,2] 3.952 0.206 3.559 3.946 4.364
beta0_pH[10,2] 3.996 0.131 3.753 3.990 4.265
beta0_pH[11,2] 1.885 0.009 1.868 1.885 1.903
beta0_pH[12,2] 1.794 0.069 1.663 1.792 1.932
beta0_pH[13,2] 2.496 0.026 2.445 2.495 2.548
beta0_pH[14,2] 1.803 0.024 1.757 1.803 1.849
beta0_pH[15,2] 2.545 0.017 2.513 2.545 2.577
beta0_pH[16,2] 2.766 0.025 2.716 2.766 2.815
beta0_pH[1,3] 1.096 0.409 0.150 1.161 1.694
beta0_pH[2,3] 2.126 0.018 2.091 2.126 2.161
beta0_pH[3,3] 2.692 0.024 2.646 2.692 2.739
beta0_pH[4,3] 2.980 0.035 2.911 2.979 3.051
beta0_pH[5,3] 3.480 2.328 0.939 2.901 9.887
beta0_pH[6,3] 1.366 0.376 0.714 1.336 2.168
beta0_pH[7,3] 0.860 0.050 0.763 0.860 0.954
beta0_pH[8,3] 0.255 0.295 -0.277 0.295 0.609
beta0_pH[9,3] -2.278 1.059 -4.852 -2.132 -0.572
beta0_pH[10,3] -0.518 0.873 -2.580 -0.395 0.843
beta0_pH[11,3] 1.128 0.507 0.089 1.159 2.038
beta0_pH[12,3] 1.003 0.661 -0.367 0.996 2.312
beta0_pH[13,3] 1.025 0.482 0.002 1.035 1.911
beta0_pH[14,3] 0.603 0.426 -0.333 0.642 1.337
beta0_pH[15,3] 0.729 0.404 -0.098 0.736 1.466
beta0_pH[16,3] 1.215 0.626 0.035 1.194 2.477
beta0_pH[11,4] 0.783 0.764 -1.487 0.914 1.902
beta0_pH[12,4] -0.146 4.209 -17.076 0.914 2.226
beta0_pH[13,4] 0.401 1.116 -3.051 0.736 1.586
beta0_pH[14,4] 0.310 0.698 -1.491 0.435 1.276
beta0_pH[15,4] 0.486 0.732 -1.237 0.684 1.532
beta0_pH[16,4] 0.761 0.963 -1.699 0.907 2.152
beta0_pH[11,5] 2.971 0.149 2.699 2.962 3.281
beta0_pH[12,5] -0.914 0.508 -1.921 -0.897 0.053
beta0_pH[13,5] 4.892 0.417 4.122 4.879 5.764
beta0_pH[14,5] 3.909 0.363 3.257 3.887 4.691
beta0_pH[15,5] 2.286 0.124 2.059 2.280 2.550
beta0_pH[16,5] 0.247 0.064 0.160 0.235 0.406
beta1_pH[1,1] 2.585 0.486 1.658 2.562 3.665
beta1_pH[2,1] 1.552 0.377 0.981 1.521 2.312
beta1_pH[3,1] 1.828 0.455 1.113 1.762 3.063
beta1_pH[4,1] 1.874 0.440 1.180 1.813 2.976
beta1_pH[5,1] 2.298 0.546 1.450 2.227 3.670
beta1_pH[6,1] 2.145 0.856 1.061 1.886 4.443
beta1_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,1] 2.327 0.933 1.167 2.073 4.667
beta1_pH[9,1] 2.061 0.670 1.171 1.929 4.019
beta1_pH[10,1] 2.149 0.468 1.460 2.067 3.300
beta1_pH[11,1] 3.372 0.548 2.517 3.303 4.733
beta1_pH[12,1] 2.897 0.279 2.363 2.887 3.466
beta1_pH[13,1] 3.452 0.245 2.988 3.449 3.945
beta1_pH[14,1] 2.840 0.338 2.202 2.826 3.572
beta1_pH[15,1] 2.397 0.308 1.900 2.369 3.013
beta1_pH[16,1] 4.037 0.305 3.462 4.024 4.639
beta1_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[1,3] 1.902 1.250 0.041 1.755 5.208
beta1_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,3] 4.659 2.401 0.920 4.160 10.655
beta1_pH[9,3] 7.584 2.660 3.586 7.144 13.822
beta1_pH[10,3] 6.617 2.726 2.823 6.048 12.825
beta1_pH[11,3] 2.157 0.551 1.082 2.160 3.263
beta1_pH[12,3] 3.217 0.681 1.947 3.186 4.644
beta1_pH[13,3] 1.953 0.519 0.951 1.949 3.023
beta1_pH[14,3] 2.168 0.486 1.300 2.128 3.225
beta1_pH[15,3] 2.096 0.474 1.195 2.084 3.070
beta1_pH[16,3] 3.660 0.685 2.216 3.680 4.916
beta1_pH[11,4] 2.010 0.792 0.715 1.918 4.221
beta1_pH[12,4] 4.233 4.216 1.802 3.160 20.983
beta1_pH[13,4] 1.593 1.118 0.313 1.290 5.004
beta1_pH[14,4] 2.056 0.769 0.905 1.942 3.976
beta1_pH[15,4] 1.515 0.765 0.323 1.359 3.261
beta1_pH[16,4] 3.619 1.013 2.104 3.460 6.231
beta1_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,1] 1.237 0.850 0.329 0.914 3.388
beta2_pH[2,1] 1.681 0.960 0.262 1.520 3.421
beta2_pH[3,1] 1.359 0.847 0.193 1.215 2.950
beta2_pH[4,1] 1.168 0.802 0.223 0.916 2.932
beta2_pH[5,1] 2.402 2.413 0.257 1.713 8.942
beta2_pH[6,1] 1.312 2.231 0.108 0.468 8.151
beta2_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,1] 1.108 1.841 0.114 0.390 7.392
beta2_pH[9,1] 1.289 1.889 0.146 0.560 7.832
beta2_pH[10,1] 1.504 1.805 0.227 0.854 7.287
beta2_pH[11,1] 0.415 0.212 0.179 0.366 0.951
beta2_pH[12,1] 1.569 0.816 0.625 1.309 3.674
beta2_pH[13,1] 1.013 0.439 0.552 0.905 2.364
beta2_pH[14,1] 1.160 0.686 0.431 0.949 3.245
beta2_pH[15,1] 1.049 0.600 0.357 0.894 2.796
beta2_pH[16,1] 1.055 0.398 0.609 0.968 1.997
beta2_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,3] 6.602 6.597 0.091 5.072 21.568
beta2_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,3] 1.557 1.462 0.168 1.075 6.155
beta2_pH[9,3] 0.301 0.605 0.087 0.200 0.988
beta2_pH[10,3] 0.437 0.752 0.101 0.267 1.967
beta2_pH[11,3] -6.262 4.915 -18.723 -4.918 -0.607
beta2_pH[12,3] -7.095 4.843 -20.189 -5.877 -1.204
beta2_pH[13,3] -6.532 4.865 -19.717 -5.244 -1.111
beta2_pH[14,3] -6.664 4.868 -19.994 -5.205 -1.256
beta2_pH[15,3] -6.770 4.816 -19.696 -5.453 -1.150
beta2_pH[16,3] -6.594 4.704 -19.623 -5.110 -1.569
beta2_pH[11,4] -5.831 4.673 -17.967 -4.368 -0.705
beta2_pH[12,4] -6.030 4.558 -17.637 -4.641 -1.047
beta2_pH[13,4] -6.010 4.502 -17.306 -4.694 -0.891
beta2_pH[14,4] -5.892 4.435 -17.086 -4.560 -0.975
beta2_pH[15,4] -5.892 4.638 -17.482 -4.416 -0.785
beta2_pH[16,4] -6.069 4.384 -17.450 -4.666 -1.260
beta2_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta3_pH[1,1] 34.964 1.166 32.822 34.932 37.571
beta3_pH[2,1] 32.921 1.410 30.446 32.839 35.931
beta3_pH[3,1] 33.053 1.619 29.984 32.946 37.087
beta3_pH[4,1] 32.790 1.944 29.509 32.737 36.917
beta3_pH[5,1] 27.625 1.151 25.931 27.493 30.731
beta3_pH[6,1] 33.278 3.203 26.471 33.142 40.548
beta3_pH[7,1] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,1] 34.309 3.526 28.742 33.897 41.704
beta3_pH[9,1] 31.955 2.817 27.153 31.625 39.698
beta3_pH[10,1] 32.752 1.507 30.207 32.648 35.872
beta3_pH[11,1] 31.412 1.234 29.292 31.368 33.954
beta3_pH[12,1] 30.783 0.471 29.792 30.796 31.636
beta3_pH[13,1] 31.424 0.527 30.434 31.398 32.594
beta3_pH[14,1] 32.500 0.753 31.236 32.445 34.144
beta3_pH[15,1] 31.028 0.716 29.534 31.071 32.348
beta3_pH[16,1] 30.122 0.375 29.421 30.117 30.854
beta3_pH[1,2] 34.505 5.500 25.488 34.340 43.496
beta3_pH[2,2] 34.337 5.461 25.395 34.284 43.494
beta3_pH[3,2] 30.610 7.308 18.521 30.649 42.396
beta3_pH[4,2] 30.438 7.191 18.575 30.352 42.386
beta3_pH[5,2] 30.483 7.276 18.554 30.700 42.317
beta3_pH[6,2] 30.526 7.255 18.609 30.672 42.400
beta3_pH[7,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[9,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[10,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[11,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[12,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[13,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[14,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[15,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[16,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[1,3] 38.067 3.494 29.280 39.471 42.696
beta3_pH[2,3] 34.586 5.458 25.509 34.742 43.530
beta3_pH[3,3] 30.436 7.245 18.629 30.300 42.398
beta3_pH[4,3] 30.460 7.207 18.617 30.630 42.465
beta3_pH[5,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[6,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[7,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,3] 42.013 1.295 38.783 42.058 43.857
beta3_pH[9,3] 36.138 1.849 34.531 35.410 41.612
beta3_pH[10,3] 37.949 2.580 34.668 37.220 43.369
beta3_pH[11,3] 40.736 1.085 38.895 40.748 42.714
beta3_pH[12,3] 42.009 0.740 40.296 42.268 42.851
beta3_pH[13,3] 41.413 0.731 39.959 41.419 42.794
beta3_pH[14,3] 40.477 0.656 39.132 40.497 41.756
beta3_pH[15,3] 40.625 0.637 39.191 40.659 41.816
beta3_pH[16,3] 41.838 0.388 41.021 41.897 42.590
beta3_pH[11,4] 40.624 1.240 38.192 40.520 43.369
beta3_pH[12,4] 41.902 0.939 40.206 42.116 43.764
beta3_pH[13,4] 42.012 1.198 39.582 42.048 43.890
beta3_pH[14,4] 40.528 0.812 38.948 40.556 42.090
beta3_pH[15,4] 41.309 1.200 38.831 41.266 43.706
beta3_pH[16,4] 41.772 0.483 40.923 41.802 42.823
beta3_pH[11,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[12,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[13,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[14,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[15,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[16,5] 25.000 0.000 25.000 25.000 25.000
beta4_pH[1,1] 0.263 0.021 0.224 0.263 0.305
beta4_pH[2,1] 0.211 0.005 0.201 0.211 0.221
beta4_pH[3,1] 0.173 0.007 0.159 0.172 0.187
beta4_pH[4,1] 0.294 0.011 0.273 0.294 0.316
beta4_pH[5,1] 0.424 0.259 0.089 0.397 0.938
beta4_pH[6,1] 0.794 0.506 0.313 0.732 1.518
beta4_pH[7,1] 0.859 0.331 0.355 0.810 1.628
beta4_pH[8,1] 0.158 0.069 0.044 0.152 0.313
beta4_pH[9,1] 0.213 0.085 0.060 0.209 0.387
beta4_pH[10,1] 0.001 0.001 0.000 0.001 0.003
beta4_pH[11,1] 0.446 0.012 0.424 0.446 0.471
beta4_pH[12,1] 0.641 0.023 0.596 0.641 0.687
beta4_pH[13,1] 0.632 0.014 0.605 0.631 0.658
beta4_pH[14,1] 0.669 0.032 0.610 0.668 0.734
beta4_pH[15,1] 0.797 0.021 0.756 0.797 0.838
beta4_pH[16,1] 1.220 0.042 1.136 1.218 1.306
beta4_pH[1,2] 0.307 0.096 0.142 0.301 0.509
beta4_pH[2,2] 0.513 0.023 0.469 0.513 0.557
beta4_pH[3,2] 0.446 0.012 0.423 0.446 0.470
beta4_pH[4,2] 0.654 0.037 0.582 0.653 0.730
beta4_pH[5,2] 0.716 0.742 0.105 0.531 2.311
beta4_pH[6,2] 0.899 0.815 0.190 0.728 2.604
beta4_pH[7,2] 0.890 0.429 0.351 0.809 1.874
beta4_pH[8,2] 0.844 0.203 0.485 0.835 1.262
beta4_pH[9,2] 1.213 1.294 0.300 0.898 4.249
beta4_pH[10,2] 0.455 0.052 0.358 0.454 0.561
beta4_pH[11,2] 1.270 0.050 1.177 1.269 1.370
beta4_pH[12,2] 0.579 0.108 0.386 0.573 0.810
beta4_pH[13,2] 0.947 0.044 0.861 0.946 1.031
beta4_pH[14,2] 0.933 0.063 0.813 0.933 1.063
beta4_pH[15,2] 1.011 0.026 0.965 1.010 1.065
beta4_pH[16,2] 1.408 0.075 1.267 1.405 1.558
beta4_pH[1,3] 0.285 0.089 0.125 0.281 0.476
beta4_pH[2,3] 0.052 0.017 0.019 0.051 0.087
beta4_pH[3,3] 0.052 0.009 0.034 0.052 0.071
beta4_pH[4,3] 0.173 0.032 0.112 0.173 0.236
beta4_pH[5,3] 0.754 0.793 0.120 0.570 2.619
beta4_pH[6,3] 0.885 0.913 0.206 0.706 2.562
beta4_pH[7,3] 0.888 0.403 0.343 0.814 1.918
beta4_pH[8,3] 0.415 0.199 0.105 0.387 0.887
beta4_pH[9,3] 0.974 1.061 0.164 0.722 3.108
beta4_pH[10,3] 0.061 0.042 0.003 0.055 0.157
beta4_pH[11,3] 0.363 0.016 0.333 0.363 0.395
beta4_pH[12,3] 0.439 0.042 0.363 0.438 0.525
beta4_pH[13,3] 0.838 0.037 0.766 0.839 0.911
beta4_pH[14,3] 0.827 0.068 0.700 0.825 0.963
beta4_pH[15,3] 0.635 0.017 0.602 0.636 0.668
beta4_pH[16,3] 0.639 0.035 0.574 0.638 0.712
beta4_pH[1,4] 0.633 0.328 0.210 0.570 1.413
beta4_pH[2,4] 0.798 0.852 0.116 0.554 3.031
beta4_pH[3,4] 0.825 0.785 0.116 0.596 2.957
beta4_pH[4,4] 0.721 0.459 0.207 0.615 1.894
beta4_pH[5,4] 0.755 0.829 0.115 0.557 2.553
beta4_pH[6,4] 0.919 0.793 0.222 0.731 2.891
beta4_pH[7,4] 0.893 0.469 0.344 0.818 1.877
beta4_pH[8,4] 0.739 0.603 0.158 0.583 2.277
beta4_pH[9,4] 0.780 0.875 0.110 0.567 2.862
beta4_pH[10,4] 1.519 3.206 0.007 0.387 11.651
beta4_pH[11,4] 0.444 0.017 0.411 0.444 0.477
beta4_pH[12,4] 0.562 0.055 0.457 0.561 0.674
beta4_pH[13,4] 0.926 0.036 0.856 0.927 0.998
beta4_pH[14,4] 0.923 0.076 0.781 0.919 1.083
beta4_pH[15,4] 0.714 0.020 0.675 0.714 0.754
beta4_pH[16,4] 0.687 0.042 0.610 0.687 0.771
beta4_pH[1,5] 0.642 0.344 0.212 0.577 1.508
beta4_pH[2,5] 0.785 0.785 0.109 0.578 2.666
beta4_pH[3,5] 0.848 0.895 0.121 0.598 3.112
beta4_pH[4,5] 0.719 0.458 0.207 0.616 1.798
beta4_pH[5,5] 0.761 0.798 0.110 0.565 2.579
beta4_pH[6,5] 0.888 0.768 0.207 0.720 2.704
beta4_pH[7,5] 0.890 0.462 0.342 0.815 1.825
beta4_pH[8,5] 0.728 0.643 0.146 0.597 2.096
beta4_pH[9,5] 0.813 1.062 0.130 0.569 2.891
beta4_pH[10,5] 1.493 3.120 0.006 0.383 11.405
beta4_pH[11,5] 0.603 0.052 0.504 0.602 0.710
beta4_pH[12,5] 0.587 0.169 0.313 0.569 0.982
beta4_pH[13,5] 0.775 0.062 0.662 0.772 0.904
beta4_pH[14,5] 0.843 0.117 0.630 0.838 1.090
beta4_pH[15,5] 1.145 0.073 1.011 1.143 1.291
beta4_pH[16,5] 18.912 3.926 10.908 19.246 24.721
beta0_pelagic[1] 2.268 0.036 2.202 2.267 2.342
beta0_pelagic[2] 1.537 0.007 1.524 1.538 1.551
beta0_pelagic[3] 0.258 0.206 -0.256 0.296 0.571
beta0_pelagic[4] 0.175 0.287 -0.577 0.216 0.648
beta0_pelagic[5] 1.443 0.155 1.157 1.438 1.752
beta0_pelagic[6] 1.353 0.059 1.240 1.353 1.468
beta0_pelagic[7] 1.479 0.039 1.405 1.478 1.559
beta0_pelagic[8] 1.672 0.043 1.590 1.672 1.757
beta0_pelagic[9] 2.796 0.065 2.673 2.794 2.934
beta0_pelagic[10] 2.603 0.031 2.545 2.603 2.665
beta0_pelagic[11] -0.137 0.188 -0.604 -0.111 0.148
beta0_pelagic[12] 1.801 0.036 1.734 1.800 1.875
beta0_pelagic[13] -1.002 0.447 -1.706 -1.040 -0.095
beta0_pelagic[14] -0.274 0.165 -0.664 -0.259 0.014
beta0_pelagic[15] -0.210 0.040 -0.323 -0.195 -0.180
beta0_pelagic[16] -0.452 0.384 -1.340 -0.361 0.071
beta1_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[3] 0.690 0.230 0.299 0.671 1.268
beta1_pelagic[4] 0.878 0.318 0.355 0.845 1.716
beta1_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[11] 1.140 0.248 0.731 1.111 1.764
beta1_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[13] 2.782 0.607 1.452 2.794 3.880
beta1_pelagic[14] 1.781 0.430 1.123 1.707 2.834
beta1_pelagic[15] 0.077 0.133 0.000 0.013 0.472
beta1_pelagic[16] 1.531 0.554 0.743 1.433 2.798
beta2_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[3] 7.045 6.279 0.280 5.413 22.823
beta2_pelagic[4] 7.268 6.306 0.232 5.863 22.633
beta2_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[11] 1.848 1.496 0.308 1.350 5.705
beta2_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[13] 0.185 0.079 0.095 0.165 0.393
beta2_pelagic[14] 0.868 0.863 0.208 0.596 3.513
beta2_pelagic[15] 2.628 1.914 0.156 2.216 7.115
beta2_pelagic[16] 0.817 1.025 0.131 0.401 4.188
beta3_pelagic[1] 28.785 3.903 22.343 28.826 35.207
beta3_pelagic[2] 28.828 3.855 22.444 28.838 35.186
beta3_pelagic[3] 28.353 2.572 23.630 29.087 32.278
beta3_pelagic[4] 25.942 1.511 22.626 26.198 28.957
beta3_pelagic[5] 29.459 6.670 18.630 29.406 40.487
beta3_pelagic[6] 29.512 6.639 18.517 29.434 40.372
beta3_pelagic[7] 29.462 6.604 18.511 29.222 40.516
beta3_pelagic[8] 29.512 6.758 18.453 29.568 40.566
beta3_pelagic[9] 29.460 6.798 18.527 29.520 40.545
beta3_pelagic[10] 29.618 6.716 18.498 29.707 40.537
beta3_pelagic[11] 30.107 1.192 27.755 30.045 32.697
beta3_pelagic[12] 30.000 0.000 30.000 30.000 30.000
beta3_pelagic[13] 31.011 4.080 24.808 30.302 40.218
beta3_pelagic[14] 34.285 1.717 31.420 34.121 38.255
beta3_pelagic[15] 33.842 5.745 23.479 34.289 41.666
beta3_pelagic[16] 31.052 2.461 26.257 31.096 35.934
beta0_yellow[1] -0.485 0.079 -0.644 -0.485 -0.325
beta0_yellow[2] 0.541 0.082 0.367 0.550 0.678
beta0_yellow[3] -0.212 0.129 -0.469 -0.217 0.014
beta0_yellow[4] 0.889 0.228 0.304 0.937 1.188
beta0_yellow[5] 0.253 0.328 -0.454 0.275 0.826
beta0_yellow[6] 0.815 0.108 0.601 0.815 1.020
beta0_yellow[7] 0.994 0.098 0.790 0.994 1.183
beta0_yellow[8] 1.026 0.094 0.843 1.027 1.210
beta0_yellow[9] 0.459 0.102 0.260 0.460 0.653
beta0_yellow[10] 0.540 0.071 0.399 0.540 0.675
beta0_yellow[11] -0.204 0.129 -0.488 -0.194 0.020
beta0_yellow[12] -1.081 0.093 -1.261 -1.081 -0.898
beta0_yellow[13] -1.356 0.206 -1.779 -1.333 -1.056
beta0_yellow[14] -0.298 0.350 -1.341 -0.236 -0.010
beta0_yellow[15] -1.253 0.164 -1.603 -1.230 -1.032
beta0_yellow[16] -0.243 0.093 -0.423 -0.244 -0.060
beta1_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[2] 1.014 0.234 0.582 1.013 1.488
beta1_yellow[3] 0.477 0.260 0.021 0.491 0.932
beta1_yellow[4] 1.329 0.603 0.643 1.174 3.155
beta1_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[11] 0.296 0.358 0.001 0.125 1.153
beta1_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[13] 1.808 1.666 0.045 1.065 6.006
beta1_yellow[14] 0.762 1.030 0.013 0.456 4.153
beta1_yellow[15] 0.170 0.232 0.011 0.076 0.796
beta1_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[2] -7.546 5.711 -21.353 -6.328 -0.398
beta2_yellow[3] -7.708 5.818 -22.184 -6.515 -0.429
beta2_yellow[4] -5.507 5.803 -20.018 -3.484 -0.127
beta2_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[11] -7.798 5.892 -22.129 -6.676 -0.361
beta2_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[13] -5.331 5.231 -18.520 -3.379 -0.275
beta2_yellow[14] -6.958 5.557 -20.018 -5.785 -0.121
beta2_yellow[15] -7.451 5.438 -20.028 -6.408 -0.298
beta2_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta3_yellow[1] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[2] 28.884 1.121 27.340 28.669 31.919
beta3_yellow[3] 32.388 3.304 23.384 32.568 38.907
beta3_yellow[4] 30.050 3.338 23.357 29.204 35.667
beta3_yellow[5] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[6] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[7] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[8] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[9] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[10] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[11] 35.065 3.671 30.314 34.078 43.333
beta3_yellow[12] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[13] 32.678 2.913 29.100 32.484 39.076
beta3_yellow[14] 33.971 3.617 29.237 33.465 42.953
beta3_yellow[15] 36.572 3.823 30.302 36.317 43.559
beta3_yellow[16] 35.000 0.000 35.000 35.000 35.000
beta0_black[1] 0.013 0.075 -0.133 0.012 0.160
beta0_black[2] 1.856 0.033 1.789 1.855 1.921
beta0_black[3] 1.237 0.037 1.166 1.237 1.311
beta0_black[4] 2.426 0.048 2.334 2.426 2.523
beta0_black[5] 5.188 2.352 1.906 4.716 11.142
beta0_black[6] 5.221 2.139 2.615 4.698 10.820
beta0_black[7] 4.095 2.188 1.612 3.509 10.131
beta0_black[8] 0.874 0.214 0.471 0.868 1.325
beta0_black[9] 2.802 0.269 2.283 2.798 3.345
beta0_black[10] 1.286 0.052 1.183 1.287 1.388
beta0_black[11] 2.979 0.051 2.878 2.979 3.078
beta0_black[12] 4.910 0.189 4.541 4.907 5.280
beta0_black[13] 0.436 0.333 0.018 0.353 1.223
beta0_black[14] 2.853 0.105 2.654 2.853 3.057
beta0_black[15] 0.723 0.053 0.620 0.723 0.825
beta0_black[16] 3.774 0.073 3.635 3.772 3.917
beta2_black[1] 0.000 0.000 0.000 0.000 0.000
beta2_black[2] 0.000 0.000 0.000 0.000 0.000
beta2_black[3] 0.000 0.000 0.000 0.000 0.000
beta2_black[4] 0.000 0.000 0.000 0.000 0.000
beta2_black[5] 0.000 0.000 0.000 0.000 0.000
beta2_black[6] 0.000 0.000 0.000 0.000 0.000
beta2_black[7] 0.000 0.000 0.000 0.000 0.000
beta2_black[8] 0.000 0.000 0.000 0.000 0.000
beta2_black[9] 0.000 0.000 0.000 0.000 0.000
beta2_black[10] 0.000 0.000 0.000 0.000 0.000
beta2_black[11] 0.000 0.000 0.000 0.000 0.000
beta2_black[12] 0.000 0.000 0.000 0.000 0.000
beta2_black[13] -2.502 1.759 -6.994 -2.104 -0.286
beta2_black[14] 0.000 0.000 0.000 0.000 0.000
beta2_black[15] 0.000 0.000 0.000 0.000 0.000
beta2_black[16] 0.000 0.000 0.000 0.000 0.000
beta3_black[1] 25.000 0.000 25.000 25.000 25.000
beta3_black[2] 25.000 0.000 25.000 25.000 25.000
beta3_black[3] 25.000 0.000 25.000 25.000 25.000
beta3_black[4] 25.000 0.000 25.000 25.000 25.000
beta3_black[5] 25.000 0.000 25.000 25.000 25.000
beta3_black[6] 25.000 0.000 25.000 25.000 25.000
beta3_black[7] 25.000 0.000 25.000 25.000 25.000
beta3_black[8] 25.000 0.000 25.000 25.000 25.000
beta3_black[9] 25.000 0.000 25.000 25.000 25.000
beta3_black[10] 25.000 0.000 25.000 25.000 25.000
beta3_black[11] 25.000 0.000 25.000 25.000 25.000
beta3_black[12] 25.000 0.000 25.000 25.000 25.000
beta3_black[13] 37.402 3.242 27.764 38.103 40.941
beta3_black[14] 25.000 0.000 25.000 25.000 25.000
beta3_black[15] 25.000 0.000 25.000 25.000 25.000
beta3_black[16] 25.000 0.000 25.000 25.000 25.000
beta4_black[1] -0.328 0.091 -0.502 -0.327 -0.152
beta4_black[2] 0.258 0.054 0.152 0.258 0.364
beta4_black[3] -0.720 0.063 -0.843 -0.720 -0.594
beta4_black[4] 0.628 0.129 0.380 0.624 0.883
beta4_black[5] 0.794 1.502 -1.152 0.503 4.439
beta4_black[6] 0.897 1.814 -0.831 0.525 4.463
beta4_black[7] 0.776 1.448 -1.058 0.497 4.249
beta4_black[8] 0.086 0.279 -0.467 0.099 0.618
beta4_black[9] 1.036 0.810 -0.088 0.878 3.084
beta4_black[10] 0.073 0.068 -0.061 0.074 0.203
beta4_black[11] -0.609 0.064 -0.729 -0.609 -0.479
beta4_black[12] 1.234 0.297 0.675 1.228 1.843
beta4_black[13] -0.999 0.077 -1.156 -0.997 -0.854
beta4_black[14] 0.220 0.183 -0.128 0.218 0.586
beta4_black[15] 0.002 0.046 -0.089 0.002 0.094
beta4_black[16] -0.418 0.116 -0.641 -0.419 -0.190
beta0_dsr[11] -0.329 0.009 -0.348 -0.329 -0.312
beta0_dsr[12] 1.409 0.044 1.323 1.409 1.493
beta0_dsr[13] 0.155 0.021 0.115 0.155 0.197
beta0_dsr[14] -0.127 0.017 -0.160 -0.128 -0.094
beta0_dsr[15] 0.205 0.014 0.178 0.206 0.232
beta0_dsr[16] 0.157 0.015 0.126 0.157 0.187
beta1_dsr[11] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[12] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[13] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[14] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[15] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[16] 0.000 0.000 0.000 0.000 0.000
beta4_dsr[11] 0.545 0.035 0.475 0.544 0.613
beta4_dsr[12] 0.802 0.232 0.353 0.802 1.261
beta4_dsr[13] -0.079 0.039 -0.157 -0.079 -0.002
beta4_dsr[14] -0.208 0.070 -0.347 -0.209 -0.072
beta4_dsr[15] 0.216 0.022 0.175 0.216 0.258
beta4_dsr[16] -0.047 0.046 -0.138 -0.046 0.043
beta0_slope[11] -2.208 0.021 -2.250 -2.208 -2.167
beta0_slope[12] -4.989 0.186 -5.371 -4.982 -4.645
beta0_slope[13] -1.501 0.033 -1.566 -1.500 -1.434
beta0_slope[14] -2.712 0.051 -2.817 -2.711 -2.616
beta0_slope[15] -1.592 0.022 -1.634 -1.592 -1.548
beta0_slope[16] -2.649 0.043 -2.734 -2.651 -2.562
beta1_slope[11] 0.000 0.000 0.000 0.000 0.000
beta1_slope[12] 0.000 0.000 0.000 0.000 0.000
beta1_slope[13] 0.000 0.000 0.000 0.000 0.000
beta1_slope[14] 0.000 0.000 0.000 0.000 0.000
beta1_slope[15] 0.000 0.000 0.000 0.000 0.000
beta1_slope[16] 0.000 0.000 0.000 0.000 0.000
beta4_slope[11] -0.336 0.075 -0.481 -0.337 -0.192
beta4_slope[12] -4.641 0.776 -6.302 -4.597 -3.223
beta4_slope[13] 0.327 0.053 0.221 0.326 0.429
beta4_slope[14] -0.160 0.152 -0.465 -0.157 0.127
beta4_slope[15] -0.361 0.033 -0.424 -0.360 -0.295
beta4_slope[16] -0.142 0.102 -0.350 -0.142 0.052
beta5_slope[11] 3.067 0.107 2.862 3.065 3.288
beta5_slope[12] 1.999 0.383 1.310 1.980 2.806
beta5_slope[13] 0.883 0.080 0.727 0.884 1.037
beta5_slope[14] 3.759 0.353 3.114 3.749 4.483
beta5_slope[15] 1.306 0.073 1.164 1.304 1.455
beta5_slope[16] 2.735 0.236 2.295 2.728 3.191
tau_comp[1] 3.766 0.300 3.201 3.754 4.403
tau_comp[2] 2.863 0.257 2.379 2.857 3.374
tau_comp[3] 0.870 0.074 0.733 0.868 1.016
tau_comp[4] 1.778 0.276 1.295 1.752 2.404
tau_comp[5] 1.262 0.189 0.917 1.253 1.665
sigma_H[1] 0.206 0.050 0.118 0.202 0.314
sigma_H[2] 0.175 0.028 0.127 0.172 0.237
sigma_H[3] 0.192 0.041 0.119 0.189 0.278
sigma_H[4] 0.461 0.080 0.328 0.452 0.644
sigma_H[5] 0.960 0.203 0.593 0.949 1.378
sigma_H[6] 0.520 0.187 0.160 0.511 0.909
sigma_H[7] 0.311 0.062 0.215 0.302 0.457
sigma_H[8] 0.417 0.087 0.274 0.408 0.608
sigma_H[9] 0.513 0.119 0.331 0.498 0.787
sigma_H[10] 0.211 0.039 0.144 0.206 0.298
sigma_H[11] 0.266 0.044 0.194 0.261 0.361
sigma_H[12] 0.382 0.147 0.203 0.333 0.724
sigma_H[13] 0.217 0.037 0.155 0.214 0.299
sigma_H[14] 0.499 0.090 0.346 0.490 0.698
sigma_H[15] 0.232 0.037 0.170 0.228 0.314
sigma_H[16] 0.216 0.041 0.150 0.212 0.309
lambda_H[1] 2.232 3.442 0.078 1.035 12.502
lambda_H[2] 5.018 4.694 0.460 3.615 17.385
lambda_H[3] 5.710 8.714 0.258 2.791 29.423
lambda_H[4] 0.006 0.004 0.001 0.005 0.018
lambda_H[5] 3.147 9.218 0.030 0.693 21.190
lambda_H[6] 6.287 14.043 0.009 0.207 46.884
lambda_H[7] 0.012 0.008 0.002 0.010 0.034
lambda_H[8] 7.949 10.190 0.118 4.414 36.763
lambda_H[9] 0.015 0.010 0.003 0.013 0.040
lambda_H[10] 0.228 0.363 0.026 0.156 0.839
lambda_H[11] 0.242 0.422 0.010 0.114 1.147
lambda_H[12] 3.674 4.879 0.167 2.096 18.229
lambda_H[13] 3.131 3.061 0.278 2.219 10.836
lambda_H[14] 3.476 4.156 0.249 2.185 14.993
lambda_H[15] 0.019 0.025 0.003 0.014 0.061
lambda_H[16] 0.466 0.748 0.031 0.246 2.331
mu_lambda_H[1] 4.068 1.902 1.000 3.871 8.133
mu_lambda_H[2] 3.674 1.957 0.507 3.467 7.813
mu_lambda_H[3] 3.382 1.852 0.683 3.116 7.600
sigma_lambda_H[1] 8.079 4.341 1.666 7.238 18.319
sigma_lambda_H[2] 8.138 4.737 0.794 7.512 18.273
sigma_lambda_H[3] 6.158 4.020 0.886 5.253 16.120
beta_H[1,1] 6.491 1.330 3.249 6.715 8.408
beta_H[2,1] 9.756 0.625 8.432 9.792 10.908
beta_H[3,1] 7.840 0.823 5.887 7.935 9.180
beta_H[4,1] 10.632 7.579 -4.885 10.613 25.119
beta_H[5,1] -0.254 2.558 -5.644 -0.081 4.079
beta_H[6,1] 2.999 4.301 -7.438 4.430 8.226
beta_H[7,1] 1.467 5.674 -10.511 1.875 11.549
beta_H[8,1] 1.058 2.117 -2.661 1.185 3.379
beta_H[9,1] 13.674 5.582 2.746 13.617 24.629
beta_H[10,1] 6.847 1.854 2.967 6.913 10.294
beta_H[11,1] 4.992 3.542 -2.901 5.741 10.003
beta_H[12,1] 2.581 1.085 0.647 2.497 4.892
beta_H[13,1] 9.096 0.897 7.214 9.153 10.653
beta_H[14,1] 2.143 1.008 0.257 2.118 4.143
beta_H[15,1] -7.051 3.380 -13.149 -7.154 0.155
beta_H[16,1] 3.983 2.979 -0.802 3.594 10.649
beta_H[1,2] 7.776 0.268 7.220 7.792 8.278
beta_H[2,2] 9.886 0.163 9.560 9.884 10.206
beta_H[3,2] 8.795 0.201 8.407 8.793 9.194
beta_H[4,2] 2.994 1.444 0.294 2.994 5.915
beta_H[5,2] 1.772 1.003 -0.145 1.764 3.739
beta_H[6,2] 6.039 1.107 3.536 6.191 7.800
beta_H[7,2] 2.377 1.073 0.494 2.306 4.623
beta_H[8,2] 3.053 0.784 1.474 3.137 4.198
beta_H[9,2] 3.153 1.104 1.000 3.133 5.399
beta_H[10,2] 8.028 0.369 7.267 8.039 8.719
beta_H[11,2] 9.814 0.638 8.857 9.705 11.273
beta_H[12,2] 3.938 0.373 3.221 3.927 4.699
beta_H[13,2] 9.162 0.245 8.695 9.155 9.653
beta_H[14,2] 4.031 0.343 3.393 4.018 4.749
beta_H[15,2] 11.556 0.612 10.287 11.565 12.663
beta_H[16,2] 4.463 0.834 2.924 4.418 6.112
beta_H[1,3] 8.444 0.264 7.965 8.435 8.982
beta_H[2,3] 9.865 0.150 9.576 9.863 10.172
beta_H[3,3] 9.455 0.167 9.145 9.451 9.793
beta_H[4,3] -2.167 0.849 -3.897 -2.172 -0.572
beta_H[5,3] 3.900 0.629 2.602 3.912 5.126
beta_H[6,3] 8.406 1.148 6.745 8.151 10.760
beta_H[7,3] -2.566 0.627 -3.870 -2.559 -1.372
beta_H[8,3] 5.266 0.419 4.675 5.203 6.307
beta_H[9,3] -2.698 0.725 -4.140 -2.709 -1.296
beta_H[10,3] 8.553 0.277 8.022 8.549 9.108
beta_H[11,3] 8.503 0.295 7.868 8.527 9.010
beta_H[12,3] 5.265 0.321 4.512 5.306 5.809
beta_H[13,3] 8.858 0.176 8.496 8.859 9.199
beta_H[14,3] 5.781 0.275 5.181 5.803 6.255
beta_H[15,3] 10.283 0.298 9.702 10.281 10.879
beta_H[16,3] 5.877 0.595 4.668 5.914 6.943
beta_H[1,4] 8.101 0.185 7.714 8.110 8.442
beta_H[2,4] 9.850 0.134 9.574 9.855 10.104
beta_H[3,4] 9.947 0.154 9.605 9.961 10.223
beta_H[4,4] 11.815 0.465 10.881 11.826 12.714
beta_H[5,4] 5.585 0.746 4.357 5.500 7.311
beta_H[6,4] 7.166 0.917 5.284 7.359 8.478
beta_H[7,4] 8.551 0.344 7.884 8.551 9.227
beta_H[8,4] 6.726 0.227 6.264 6.732 7.140
beta_H[9,4] 7.429 0.449 6.560 7.418 8.336
beta_H[10,4] 7.526 0.224 7.099 7.518 7.982
beta_H[11,4] 9.511 0.192 9.123 9.516 9.888
beta_H[12,4] 7.248 0.192 6.881 7.249 7.623
beta_H[13,4] 9.067 0.140 8.784 9.071 9.338
beta_H[14,4] 7.834 0.218 7.428 7.831 8.272
beta_H[15,4] 9.473 0.221 9.034 9.473 9.899
beta_H[16,4] 9.435 0.239 8.977 9.429 9.912
beta_H[1,5] 8.894 0.158 8.572 8.901 9.186
beta_H[2,5] 10.679 0.101 10.493 10.675 10.896
beta_H[3,5] 10.755 0.164 10.456 10.748 11.086
beta_H[4,5] 8.347 0.490 7.411 8.326 9.335
beta_H[5,5] 5.189 0.597 3.699 5.262 6.151
beta_H[6,5] 8.981 0.614 8.051 8.870 10.308
beta_H[7,5] 6.577 0.341 5.910 6.569 7.273
beta_H[8,5] 8.241 0.193 7.917 8.223 8.643
beta_H[9,5] 8.187 0.462 7.249 8.186 9.111
beta_H[10,5] 10.047 0.230 9.579 10.048 10.490
beta_H[11,5] 11.594 0.221 11.162 11.591 12.036
beta_H[12,5] 8.610 0.190 8.240 8.606 8.990
beta_H[13,5] 10.142 0.130 9.888 10.138 10.414
beta_H[14,5] 9.234 0.216 8.820 9.227 9.695
beta_H[15,5] 11.315 0.233 10.858 11.315 11.769
beta_H[16,5] 9.828 0.186 9.437 9.835 10.162
beta_H[1,6] 10.123 0.244 9.756 10.085 10.715
beta_H[2,6] 11.468 0.117 11.235 11.468 11.692
beta_H[3,6] 10.620 0.155 10.287 10.629 10.899
beta_H[4,6] 12.923 0.867 11.191 12.936 14.588
beta_H[5,6] 5.669 0.627 4.625 5.622 6.982
beta_H[6,6] 8.720 0.576 7.419 8.752 9.781
beta_H[7,6] 10.566 0.589 9.351 10.570 11.702
beta_H[8,6] 9.534 0.246 9.061 9.539 9.967
beta_H[9,6] 8.289 0.768 6.795 8.263 9.850
beta_H[10,6] 9.136 0.340 8.429 9.159 9.762
beta_H[11,6] 10.688 0.359 9.917 10.720 11.308
beta_H[12,6] 9.368 0.241 8.887 9.360 9.871
beta_H[13,6] 11.094 0.154 10.801 11.091 11.412
beta_H[14,6] 9.921 0.286 9.308 9.932 10.453
beta_H[15,6] 10.550 0.417 9.713 10.548 11.344
beta_H[16,6] 10.583 0.260 10.040 10.594 11.066
beta_H[1,7] 10.134 1.265 6.780 10.423 11.883
beta_H[2,7] 12.123 0.521 10.980 12.145 13.090
beta_H[3,7] 10.299 0.671 8.783 10.355 11.394
beta_H[4,7] 1.401 4.503 -7.282 1.372 10.710
beta_H[5,7] 5.850 2.049 2.148 5.899 10.031
beta_H[6,7] 7.052 3.178 -0.977 8.150 10.735
beta_H[7,7] 7.597 2.995 1.751 7.562 13.794
beta_H[8,7] 10.854 0.838 9.293 10.861 12.350
beta_H[9,7] 4.959 3.930 -3.126 5.097 12.647
beta_H[10,7] 10.057 1.568 7.252 9.934 13.456
beta_H[11,7] 11.345 1.747 8.314 11.190 15.288
beta_H[12,7] 9.782 0.954 7.724 9.870 11.381
beta_H[13,7] 11.676 0.708 10.086 11.761 12.828
beta_H[14,7] 10.633 0.896 8.807 10.674 12.271
beta_H[15,7] 14.004 2.201 9.819 14.025 18.470
beta_H[16,7] 12.911 1.403 10.484 12.806 16.084
beta0_H[1] 8.948 18.927 -28.629 8.875 50.028
beta0_H[2] 10.679 8.115 -5.825 10.708 27.227
beta0_H[3] 9.873 10.128 -10.216 9.958 30.575
beta0_H[4] 10.122 180.032 -353.269 9.747 382.268
beta0_H[5] 3.774 27.192 -47.567 4.353 50.618
beta0_H[6] 7.946 51.207 -103.158 7.583 125.900
beta0_H[7] 3.632 135.887 -266.884 2.863 275.821
beta0_H[8] 6.446 21.694 -16.589 6.445 27.981
beta0_H[9] 11.208 124.940 -237.541 9.693 263.576
beta0_H[10] 8.316 38.449 -73.297 8.837 83.369
beta0_H[11] 8.201 50.439 -99.031 8.776 119.188
beta0_H[12] 6.624 12.174 -17.587 6.492 31.183
beta0_H[13] 10.427 11.635 -11.569 10.116 32.925
beta0_H[14] 7.204 11.170 -15.766 7.250 30.174
beta0_H[15] 6.138 112.883 -218.945 4.213 239.389
beta0_H[16] 7.333 33.649 -62.087 7.973 71.790